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Abstract 

We discuss the topological CP^ model which consists of the holo- 
morphic maps from Riemann surfaces onto CP^. We construct a 
large-iV matrix model which reproduces precisely the partition func- 
tion of the CP^ model at all genera of Riemann surfaces. The action 
of our matrix model has the form TrF(M) = — 2TrM(logM — 1) + 
2Etn,pTrM"(logM-c„)+El/n-tn-i,QTVM- (c„ = Ei Vj) where 
M is an A^ X A^ hermitian matrix and t^.p itn,Q), {n = 0,1,2- ■ ■) are 
the coupling constants of the n-th descendant of the puncture (Kahler) 
operator. 



The topological CP^ model was introduced a few years ago by |1], 
as a characteristic example of a topological a model coupled to the two- 
dimensional gravity. As is well-known, the topological cr-model (or topo- 
logical A model) is a twisted version of the supersymmetric non-linear cr- 
model with a target space which has an almost complex structure Q. Its 
path-integral is given by a sum over contributions from instantons which are 
holomorphic maps from the (genus g) Riemann surface onto a target space 
K. When the target space K is CP^, the theory possesses two observables 
P and Q corresponding to the identity and the Kahler class of CP^ and 
thus the theory becomes a two-primary model. When the system is coupled 
to the two-dimensional gravity, additional excitations, e.g. gravitational de- 
scendants an{P), (TniQ) {n = 1,2, ■ ■ ■) appear in the theory (by convention 
(Tq[P) = P,ao{Q) = Q). We denote the coupling constants associated with 
these observables as t^^p, tn^q {n = 0, 1, 2, ■ ■ ■). 

In the following we construct a matrix model which reproduces exactly 
the partition function of the CP^ model in the large phase space {tn„p 7^ 
0, tn,Q 7^ 0) at all genera of the Riemann surface. The potential of our 
matrix model has the form 

V{M) = -2M(logM - 1) + 2 ^ t„,pM"(logM - c„) 

ra=l 
1 " 1 

+ E -tn-l,QM- , C^ = T.- (1) 

n=l ^ j=l 3 

and the partition function is given by 

Z= f dMexp{NTTV{M)) . (2) 

Here M is an A^ x A^ hermitian matrix. The characteristic feature of our 
model (1) is the appearance of the logarithmic potentials indicative of the 
asymptotic freedom and the non-vanishing j3 function. In the limit when all 
coupling constants vanish our action is reduced to V{M) = — 2M(log M — 1) 
which is the form of the effective potential of the CP^ model when M is 
replaced by the chiral superfield $ 0]. The extremum of the action V'{^) = 
occurs at $^ = 1 which is the basic relation of the CP^ quantum cohomology 

Let us first describe some basic features of the CP^ model at genus g = 0. 
By denoting the free-energy of the system as F basic two-point functions of 



the theory are defined as 

{PP) = Ji^^n, (3) 

d'^F 
dto^qdio^p 

A crucial ingredient of the CP^ model is the following (constitutive) re- 
lation 

{QQ) = e^PP) = e" (5) 

which conies from the instanton analysis |I| . (|^) is the formula which distin- 
guishes the CP^ model from other two-primary models like the Ising model. 
From (H),(§),® we obtain the fiow equations for the time variable to,Q 

--V , (6) 



dt, 



0,1 



3v 

= (eT- (7) 



Here ' denotes the derivative in to.p- The fiow equations for the descen- 
dant times tn,p, tn,Q {n = 1,2- ■ ■) are obtained by using the genus g = 
topological recursion relations ^ 

(a„($,)XF) = n{an-i{^a)P){QXY)+n{an-i{^a)Q){PXY) , (8) 

where $q,=P or Q. du/dti^p is, for instance, derived as 

du 



dti 



Similarly 



dv 

dtZf 



{(Ji{P)PP) = {PP){QPP) + {PQ){PPP) 

uv' + vu' = {uv)' . (9) 

{a,{P)PQ) = {PP){QPQ) + {PQ){PPQ) 

M(e")' + W = (it;2 + (M-l)e")' . (10) 



For the sake of illustration we also present the flow equations for the variables 
ti,Q, t2,p and t2,Q 

^" ;i..^ + e-)', ^ = (-7, (11) 



^"^ H2 + 2(M-2)e"y, ^ = (it;3 + 2i;(m - l)e")' , (12) 



dt2,p dt2,p 3 

^"^ :it;=^ + 2t;e")' , -^ = (^;V + e^")' . (13) 



^^2,Q 3 9t2,Q 

It is clear that one can write down flow equations for any of the variables 
in,p^ tn,Q by using the recursion relations repeatedly. By construction these 
flows mutually commute with each other and thus they define an integrable 
hierarchy. 

It is easy to see that the CP^ system is closely related to the Toda lattice 
hierarchy p. In fact by combing (^ and (^ we obtain 

^^"^ =(e")" (14) 



which may be regarded as the continuum (genus g = 0) limit of the Toda 
lattice equation 

-^ = e""+i + e""-^ - 2e"" (15) 

ozoz 

if we identify to,p = n/N {N is the size of the lattice). In (0) to,Q ^ z = z 
and thus we have the one- dimensional Toda theory. As we shall see in the 
following, our system can in fact be described by a suitable generalization of 
the one-dimensional Toda lattice hierarchy. 

Let us now introduce the Lax formalism for the Toda lattice (at genus 
g = 0) and see how one can reproduce the evolution equations of the CP^ 
model. We consider an operator 

L = p + v + e>-i , (16) 

where the momentum p is the canonical conjugate of to,p- The Poisson 
bracket of the Toda system is defined as [0] 

,, „. /M dB dB dA X ,^^. 



Time-evolution is generated by the Hamiltonians {L^)+, n = 1, 2 ■ ■ -, where 
+ means taking terms with non-negative powers of p. By introducing the 
To da-times t„ we have 

^ = {{L-U,L}, n = l,2,---. (18) 

By comparing (|l^) with the flow equations of the CP^ model we find that 
the CP^-times tn^qin = 0, 1, 2 ■ ■ ■) are nothing but the Toda-times tn+i{n = 
0,1,2---) 

tn,Q = (n + l)tn+i , n = 0,1,2,- ■■ . (19) 

For instance, 

dL 

— = {L+,L} = {p + v,p + v + ey-^} 

= (e")' + v'ey-' . (20) 

Hence dv/dti = (e^)', du/dti = v' which coincide with eqs.(|^,(^). Similarly 

dL 



^^ {(L')+, L] = {p' + 2vp + v' + 2e",p + v + eV } 

012 



= 2(t;e")' + {v' + 2e")'eV . (21) 

Hence dv/dt2 = 2(fe")', du/dt2 = 2(f^/2 + e")' which agree with eq. (|TT|) . 
Thus the generators of the t„,Q-time evolutions are identified as simple powers 
of the Lax operator L. 

If one attempts to describe the t„^p-time evolutions using the Lax formal- 
ism, one has to generalize the standard treatment of the Toda theory and 
introduce the logarithm of the Lax operator logL. We introduce a series of 
the operators L"(logL — c„) with Cn = Yl]=i 1/j and consider an additional 
set of evolution equations 

dL 

— = {(L"(logL-c„))_,,L}, n = 0,l,2,---. (22) 

It is possible to show that only terms proportional to p^ and p~^ appear in 
the right-hand-side of eq. (|22| ) as in the case of (p!8| ) (for the treatment of the 
logarithm see below, eq.(p4D) and hence (^) gives rise to a set of differential 



equations for the coefficient functions u and v. Furthermore, the flows in the 
variables s„ commute among themselves and also with those of the variables 

Then by comparing p^) with the CP^ t„.p-flow equations we find 

1 . , 

2'Sn = tn,p , n = 0,1,2, •••. (23) 

Let us check the n = case. We ffist rewrite the logL operator as 

log L = log{p + V + e^p^-^) 

1 1 

= - logp(l + vp^^ + e>~^) + - logp~^e''(l + ?;e"> + e"V) 

= -u + - log(l + vp'^ + eV^) + - log(l + ve-y + e"V)- (24) 

Taking the positive part simply drops the 2nd term in (0). Then by direct 
calculation, 

^ = {(logL)^,L} = -(.;'+ (e7p-^) ■ (25) 

Hence 

So = 2to,p . (26) 

We can check other cases n = 1, 2, ■ ■ ■ of ( ^3]) with somewhat lengthy cal- 



culations. The coefficients c„ are determined to eliminate the d/dtn-i^q 
component in the action of the operator (L"'logL)_|_ on L. Our operators 
L"(logL — Cn) have a special "scaling" property 

^L"(logL - Cn) = nL^-\\ogL- c„_i) (27) 

which plays an important role in what follows. We summarize our identifi- 
cation of the CP^ evolution equations at the genus g = Q level, 

BT 

nj- = {(L")+,L}, n = l,2,---, (28) 

1 dr 

-^— = {{L"i\ogL-c^))^,L}, n = l,2,--- . (29) 

Before we turn to the subject of the matrix-model realization of the hier- 
archy (p8D,(p9D we must discuss its higher-genus corrections. As we discussed 
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elsewhere |p, it is straightforward to determine corrections to the flow equa- 
tions of the CP^ model at each genus simply by demanding the commuta- 
tivity of the flows. Up to the order of 5^ = 2 the corrected flow equations (for 
to,Q, ti,p and ti^q) are given by 



du 
dto,Q 
dv / „ , 9 . 1 ,9 1 



du 



v' , (30) 

+A^(— u"" + —u'u'" + -^u'^u" + -^u'^ + —u"^)e^)' , 
^360 180 1440 1920 240 ' ^ 

(31) 

^ = (^-v' + {u-l)e^ + \\^^{u + ?,y' + ^^{u + 2)u")e- 

+A^(^(tx + A)u"" + ^^{u + h)u'u"' + Y^(« + Q)u"u'' 

+ ^(u + 7)u'^ + — (M + 5)u"2)e"y , (33) 

1920^ ' 240^ ^ ^ y ' ^ ^ 



V9 ' '94^ ^ 



dti^Q V2 24 

+A^(— u"" + —u'u'" + — u'^m" + — m'4 + — M"2)e"y , (34) 
420 48 48 384 72 ' ^ ' ^ ' 

-^ = Le'' + ^(u'^v + 2u'v' + 2u"v + Av")e'' 
(jti,Q 24 

+\\—vu"" + — t;M'M"' + -^vu'^u" + ^t;^'^ 
^360 180 1440 1920 

+ —vu"^ + —v'u'" + — ^;'m'm" + — t;'u'3 
240 180 720 480 

+ —v"u" + — t;"u'2 + — t;'"u' + —v"")e'')' . (35) 

720 120 80 120 ^ ^ ^ ^ 

Here A^ is the genus expansion parameter. In deriving the above equations 



we first write down candidate correction terms in the right-hand-side of (pTI)- 
( |35D with unknown coefficients and then fix them by demanding the com- 
mutativity of the flows, e.g. the vanishing of the cross- derivatives such as 
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d / dti^p{dv / dto^o) — d / dto^Q^dv / dti^p) = 0. The flow commutativity leads 
to an over-determined system for the unknown coefficients and gives rise to 
their unique solution. 

We now turn to the construction of a matrix model which reproduces the 
CP^ model. Let us ffist recall some basic formulas for the analysis of matrix 
models using orthogonal polynomials 0. We consider a one-matrix model 
with an action Tr V{M) which depends only on the eigenvalues of the N x N 
hermitian matrix M. The orthogonal polynomials are defined by 

rfAe^^("Vn(A)Vn^(A) = Sr^^mhm, (36) 

^^(A) = A" + lower terms . (37) 

The multiplication by A of the orthogonal polynomials leads to a three-term 
recursion relation 

XtpnW = Yl Qnmi'mW = ^n+l(A) + t;„V^n(A) + -R„V^„_l(A) , (38) 
m 

Rn = hJK^, = e^(<^"-^"-i). (39) 

The free-energy of the system is given by F = \ogU.Q~^hn = NY.o~^ (l>n- 
The matrix Q has non-vanishing elements only on the diagonal and 1st off- 
diagonal lines and thus is a Jacobi matrix. In the continuum limit A^(0n ~ 



F" /N"^ = u and the matrix Q turns into the Lax operator 



Qr:^L=p + v + ey-^ . (40) 

Q replaces L in the discussion of the CP^ model at higher genera. Time 
evolution equations are now given by the commutation relations 

n dQ 



1 dQ 

2N dUp 



■■[iQn+,Q] , n=l,2,--- 


(41) 


[(Q"(logQ-c„))+,Q], n = 0,l,2,---, 


(42) 



where -|- means to take the upper-triangular part of a matrix including the 
diagonal line. It is easy to check that the matrix commutators reduce to the 
Poisson bracket (17) in the continuum limit and thus ( ^l]) ( ^2|) reduce to (^) 
(E9|) at genus g = 0. 



The derivative operator d/dX, on the other hand, becomes a lower- 
triangular matrix when acting on the orthogonal polynomials 

— V'nlA) = mjJn^iiX) + --.=J2 Pnm^m{X) ■ (43) 

By partial integration we have 

rfAe^'^W^7An(A)^„(A) = -|dAe^^WiVV^'(A)^„(A)^„(A) 

= -NV'{Q)nmhm , n>m. (44) 

Hence 

P = -NV'iQ)^ , (45) 

where — means to take the lower-triangular part of the matrix. The string 
equation then reads 

[Q,P] = -N[Q,V'iQ)-] = l- (46) 

Let us now introduce our matrix model defined by an action 
TrV(M) = -2TrM(logM-l) + ^2t„,pTrM'"(logM-c„) 

n=l 

+ 5]-t„_i,QTrM". (47) 



n=l'^ 



We then obtain 



N[V'{Q)^,Q] 

-2iv[(iogg)_,Q] + iv5^2nt„,p[(g'^-^(iogg-c„_i))_,g] 
+iv^t„_i,Q[(g"-i)_,g] 

- E ^tn,Pjr, E(^ - ^)tn-l,QTr^ , (48) 



where we used (PD,(^2D and the "scaling" property of our logarithmic po- 
tentials (|27|). (46) is now compactly expressed as 

l + EE-^n,.-^^ = |^, ia = P,Q). (49) 



The diagonal element of P9| ) reads 

^^ = 1 + EE^wK-i($.)pg) (50) 

in terms of the correlation functions. By integrating (|50D we obtain 

« = *0,Q + EE nt„,aK-l($a)i^) , (51) 

n a 

V = to,P + EE ntn,a{(rn-li'^a)Q) ■ (52) 

n a 

These are the string equations of the CP^ model 0] . 

Therefore, in order to demonstrate the realization of the CP^ model by 
our matrix integral we finally have to examine our flow equations (p]),(^) 
and check if they agree with the flow equations of the CP^ model at higher 
genera. Checking (^T]) is straightforward. In the cases n = 1,2, for instance, 
(^rp gives 



dto,Q 
dv„. 



dt 



Vn , (53) 

iV(K+i - Rn) , (54) 



0,Q 
^ = liRn^l + Rn + vl), (55) 

■^r^ = —ii'^n+l + Vn)Rn+l " {Vn + f„-l)i?n) • (56) 



,Q 



If we use the definition to,p = n/N,Un = 4>n,Rn = exp (A^(0„ — 0n-i)) 
and Taylor-expand Rn±i{vn±i) around -R„(f„), we easily recover the CP^ 
equations. For instance, by expanding the right-hand-side of (Kl) as 



iV(i?„+, - /2„) = e-(< + -^ + 



U„ U^ '^n'^n 



12iV2 24Ar2 e]Y2 
1920iV4 UAN^ 12N^ 96A^4 
72iV4 120iV4 360iV4 ^ ' '^ ^ '' 
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we reproduce (0). Note that the genus expansion parameter A^ becomes 



Checking (021), on the other hand, is rather non-trivial. We first have to 



compute the matrix elements of the logarithm of a Jacobi matrix Q. We use 
a formula, 

(logg)„^ = ^.i — ^"-'"log (z + z'^R{n + z^) + v{n + z^)) ■ 1 (58) 
2m J z ^ dz dz ' 

and expand the arguments of R and v. We can then compute logQ pertur- 
batively for each fixed order of differentiations of i?, v. Up to the 2nd order 
in the derivatives, for instance, we obtain 

(logQW. = liogfl„ + jL|i + 

1 1 1 

(logQ)„,„_i = i:Vn- —-v'^^ 



2 4A^ " 24A^2 



1 K 

24iV2 R^ 


1 K 

2Am Rl ' 


(59) 


<, 




(60) 



^iogc^K„+i - 2/2„ + 4iVVi?„ Rl )^2AN^Rn 

_j_Ryn j_r;^ _ J_K< (... 

4N^ Rl 2N^ Ri m^ Rl ' ^ ' 

Using the explicit expressions for log Q as above one can evaluate the right- 
hand-side of (|4^ ) and check its agreement with the CP^ equations. We have 
calculated logQ up to the 3rd order in the derivatives and have confirmed 
the agreement of (^) with the tn,p flow equations up to ra = 2 and g = 1. 

We should note that in our construction we take the ordinary large- A^ 
limit of the matrix model rather than the double scaling limit: we do not 
adjust the coupling constants to any specific set of values. This is in accord 
with the idea that the CP"^ model, due to its asymptotic freedom, has no 
phase transition points. In this sense our model is in a spirit closer to the 



Kontsevich [T^ or Penner |Tl| model rather than the double-scaled matrix 
models describing the c < 1 minimal matter coupled to the two-dimensional 
gravity [jl2[- We think it is fortunate that the CP^ model has a relatively 



simple representation in terms of a single matrix integral. Our model may 
prove useful in the study of the moduli space of holomorphic maps from 
Riemann surface to CP"^. It will be very interesting to see if our construction 
may be generalized to CP" or Grassmann manifolds. 
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